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ABSTRACT 

Aims. We investigate the frequency-dependent radio properties of the jet of the luminous high-redshift (z = 3.2) radio quasar 
PKS 1402+044 (J1405+0415) by means of radio interferometric observations. 

Methods. The observational data were obtained with the VLBI Space Observatory Programme (VSOP) at 1.6 and 5 GHz, supple- 
mented by other multi-frequency observations with the Very Long Baseline Array (VLBA; 2.3, 8.4, and 15 GHz) and the Very Large 
Array (VLA; 1.4, 5, 15, and 43 GHz). The observations span a period of 7 years. 

Results. We find that the luminous high-redshift quasar PKS 1402+044 has a pronounced "core-jet" morphology from the parsec to 
the kilo-parsec scales. The jet shows a steeper spectral index and lower brightness temperature with increasing distance from the jet 
core. The variation of brightness temperature agrees well with the shock-in-jet model. Assuming that the jet is collimated by the ambi- 
ent magnetic field, we estimate the mass of the central object as ~ 10 9 M o . The upper limit of the jet proper motion of PKS 1402+044 
is 0.03 mas yr" 1 (~ 3c) in the east-west direction. 
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1. Introduction 

Very Long Baseline Interferometry (VLBI) studies of high- 
redshift quasars at a given observing frequency v b s can fa- 
cilitate comparison of their structural properties with those of 
their lower-redshift counterparts at a higher frequency, v b s = 
v em /(l +Z), where v em is the emitted (rest-frame) frequency and z 
the redshift of a distant quasar. High-redsfhit quasars provide in- 
dispensable input in all kinds of studies of the redshift-dependent 
properties of extragalactic objects, such as the apparent "angular 
size - redshift" ("8 - z", e.g. Gurvits et al. |1999l l and "proper 
motion - redshift" ("/i - z", e.g. Kellermann et al. 1999) rela- 
tions. 

A statistical study of 151 quasars imaged with VLBI at 
5 GHz (Frey et al. 1997) demonstrated an overall trend of a de- 
creasing jet-to-core flux density ratio with increasing redshift, 
which could be explained by the difference in spectral indices of 
cores and jets. Furthermore, a majority of radio QSOs at z > 4 
seemed to be even more compact than expected from the direct 
extrapolation of the properties of quasars at z < 4 (Paragi et al. 
[19991 . 

A number of morphological studies of high-redshift objects 
have been made with the Japanese-led Space VLBI mission 
VSOP (VLBI Space Observatory Programme). Observations 
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with the VSOP utilised an array consisting of a group of Earth- 
based radio telescopes and an 8-m space-borne antenna on board 
the satellite HALCA (Hirabayashi et al. 11998b . The orbiting an- 
tenna with an apogee of ~ 21 000 km and perigee of ~ 560 km 
provided milli-arcsecond (mas) and sub-mas resolution at the 
frequencies of 1.6 and 5 GHz. VSOP observations at 1.6 GHz 
provided roughly the same angular resolution as Earth-based 
VLBI observations at 5 GHz. Thus, dual-frequency VSOP ob- 
servations made it possible to map the distribution of spectral 
index across the source structure (e.g. Lobanov et al. 20063 and 
to study frequency-dependent structural properties. 

PKS 1402+044 (J1405+0415) is a flat-spectrum radio source 
from the Parkes 2.7-GHz Survey. In optics, it is a 19.6- 
magnitude (g filter) stellar object at a redshift of z = 3.215. It is 
a weak X-ray source with count rates of (5.6 + 1.2) x 10~ 3 ct s , 
over the band 0.2 - 4 keV in the Einstein IPC survey database 
(Thompson et al. [19981 and (1.3 + 0.2) x 10~ 2 ct s _1 over 
the band 0.1 - 2.4 keV in the ROSAT observation (Siebert 
et al. 119981 1. The Multi-Element Radio Linked Interferometer 
Network (MERLIN) observations of PKS 1402+044 made at 
1 .6 GHz indicates that there is a secondary component at a sep- 
aration of 078 at a position angle of -123° and a faint extended 
emission at a distance of 3'.'3 at a position angle of -106°. VLBI 
observations at 5 GHz (Gurvits et al.Q992) found that the main 
component consists of a compact core and a resolved jet extend- 
ing up to ~ 18 mas to the west. 
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The quasar PKS 1402+044 represents a relatively rare case 
of a strong radio source at z > 3 and therefore a potentially re- 
warding target for structural studies within a broad range of an- 
gular scales. VSOP observations with their record-high angular 
resolution at 1 .6 and 5 GHz facilitate direct comparison of struc- 
tural properties of PKS 1402+044 with its more abundant strong 
radio quasars at lower redshifts at the same emitting frequency. 

In this paper, we present VSOP images at 1.6 and 5 GHz, 
a 15-GHz Very Long Baseline Array (VLBA) image, and Very 
Large Array (VLA) images at 1.4, 5, 15, and 43 GHz of the 
quasar PKS 1402+044; discuss its spectral properties and bright- 
ness temperature variation along the jet; and determine some 
physical parameters of the core and the jet. Throughout the 
paper, we define the spectral index a as S v <x v° and adopt 
the ACDM cosmological model (Riess et al. 120041 ) with H = 
75 km s" 1 Mpc~', Q m = 0.3, and Q A = 0.7. In the latter model, 
the linear scale factor for PKS 1402+044 is ~ 7 pc mas~' . 



2. Observations and data reduction 

2.1. VSOP experiment 

Using the space-borne radio telescope HALCA and the VLBA, 
we observed the quasar PKS 1402+044 in left hand circular 
polarization at 5 GHz on 20 January 2001 and at 1.6 GHz on 
21 January 2001. The observations lasted for 8 h at 1.6 GHz 
and 7 h at 5 GHz. The data were recorded using the VLBA 
tape system with a 32-MHz bandwidth consisting of 2 interme- 
diate frequency (IF) bands and 2-bit sampling, corresponding to 
the data rate of 128 Mbps. Four tracking stations (Goldstone, 
Robledo, Tidbinbilla, Green Bank) were used to receive the 
HALCA downlink data. The HALCA data were recorded for 
~5 h at 1.6 GHz and 4.2 h at 5 GHz. The data were correlated at 
the VLBA correlator in Socorro with 128 spectral channels and 
an integration time of 4.2 s for the ground-ground baselines, and 
2.1 s at 1.6 GHz, 1.0 s at 5 GHz for the space-ground baselines. 
In the 1.6-GHz observation, the Tidbinbilla station lost 55 min- 
utes of space data and Green Bank lost all the space data (34 min- 
utes) due to a problem with recording. At 5 GHz, the IF 1 data 
were lost for 40 minutes due to a technical malfunction at the 
Tidbinbilla station. Except for these problems, fringes were suc- 
cessfully detected on all the space-ground baselines at all times. 

A priori calibration was applied for both datasets using the 
Astronomical Image Processing System (AIPS; Cotton 1995). 
After correcting the amplitudes in cross-correlation spectrum us- 
ing measurements of auto-correlation spectrum and dispersive 
delay due to the ionosphere from the maps of total electron con- 
tent, we applied a priori amplitude calibration from the antenna 
gain and system temperature measurements at the Earth-based 
telescopes. We used the respective nominal values^ for HALCA. 
After inspecting the IF bandpass, the side channels (1-5, 105 - 
128) in each IF were deleted because of the lower amplitude 
(< 80%) than in the centre channels. This reduced the useful ob- 
serving bandwidth to 22.75 MHz. Some channels affected by 
radio frequency interference were flagged, too. We corrected 
the residual delays and rates using a two-step fringe-fitting. We 
first fringe-fitted the ground-array data with a solution interval 
of 2 minutes. Then we applied the solutions to the data, fixed 
the calibration for ground antennas and determined the calibra- 
tion solutions of the space antenna using fringe-fitting with a 
4-minute interval. After that, we combined all fringe solutions, 
applied them to the data, averaged all the channels in each IF, and 
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Fig. 1. The effective (u,v) coverage of the VSOP observations of 
PKS 1402+044 at 1.6 GHz (bottom) and 5 GHz (top). At each 
frequency, the inner tracks correspond to the ground-ground 
baselines, and the outer tracks denote the space-ground base- 
lines. The rectangle in the 5-GHz (u, v) coverage shows the size 
of the 1.6-GHz (u, v) coverage. The space-ground baselines pro- 
vide Mv-ranges roughly 2.5 - 3 times longer than ground-ground 
baselines. 
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split the multi-source data into single-source data sets. Finally, 
the data were exported into Difmap (Shepherd et al. 1994) and 
averaged further over 60 s time intervals. The hybrid imaging 
and self-calibration were done in Difmap. The resulting effective 
(u, v) coverages of the VSOP observations are shown in Fig.Q] 
The correlated flux densities as a function of the projected base- 
line length are displayed in Fig. [2] top and middle. 
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2.2. VLBA and VLA data 

The 15-GHz VLBA data presented in this paper are from the 
VLBA-VSOP support survey by Gurvits et al. (in preparation). 
The observations were conducted on 5 December 1998 with left- 
hand circular polarisation, 64-MHz bandwidth and ~50-minute 
on-source observing time. We also used the 2.3/8.6 GHz visi- 
bility data provided by US Naval Observatory (USNO) Radio 
Reference Frame Image Database (RRFIDjJ Those observa- 
tions used 10 VLBA antennas and some additional geodetic an- 
tennas. All the VLA data used in this paper were obtained from 
the NRAO Data Archival The basic parameters of the VLA ob- 
servations are summarised in Table Q] The columns give (1) fre- 
quency in GHz, (2) program ID, (3) date as dd/mm/yy, (4) array 
configuration, (5) antenna numbers, (6) bandwidth in MHz, and 
(7) total on-source time in seconds. All the VLA observations 
used 3C 286 as the prime flux density calibrator. After a pri- 
ori calibrations in AIPS, we performed self-calibration, imaging 
and model-fitting in Difmap. 



Table 1. VLA observations summary. 



Vobs 


Program 


Date 


Conf. 




BW 


TOS 


(GHz) 




dd/mm/yy 






(MHz) 


(s) 


1.4 


AH0633 


11/03/98 


A 


23 


100 


130 


4.8 


AG0670 


09/10/04 


A 


26 


100 


2010 


15.9 


AH0633 


11/03/98 


A 


27 


100 


170 


43.3 


AL0618 


26/01/04 


BC 


26 


100 


1320 



3. Results 

Figure [3] shows all the images of PKS 1402+044 of the cur- 
rent study. Their parameters are summarised in Table [2] The 
columns give: (1) frequency in GHz, (2) array and configura- 
tion, (3) weighting (NW: natural, UW: uniform), (4-5) size of the 
synthesised beam in mas, (6) position angle of the major axis in 
mas, (7) peak brightness in mJy/beam, and (8) image noise level 
in mJy/beam (1 cr). From the final fringe-fitted VSOP data set 
we produced two images: (1) an image with all data included 
(hereafter VSOP image) and (2) an image using only the ground 
VLBA data at each frequency. In the imaging process, we scaled 
the gridding weights inversely with the amplitude errors. The 
VSOP image fidelity is limited, in particular, by the complete- 
ness of the (u, v) coverage. In our case, the latter is essentially 
one-dimensional (see Fig.[T) that leads to a highly elliptical syn- 
thesised beam. However, luckily, the highest angular resolution 

is achieved along the position angle of 60°, very close to that 

of the inner pc-scale jet. Thus, the space-ground baselines ob- 
tained play an important role in imaging the inner pc-scale jet of 
PKS 1402+044. 

We detect a very weak jet emission extending up to ~ 
150 mas (~1 kpc) in the high dynamic range (~ 10000) VLBA 
image at 1.6 GHz show Fig. [3^. It represents a typical core-jet 
morphology. We identify a compact core (component A) and five 
jet emission regions (components B - F) using circular Gaussian 
model-fitting in Difmap. The parameters of the models are listed 
in Table [3] The columns give: (1) component identification, (2) 
total flux density of the component, (3 - 4) radius and position 
angle of the centre of the component, (5) size of the fitted cir- 
cular Gaussian model, (6) the smallest detectable size, and (7) 

2 http://rorf.usno.navy.mil/RRFID 

3 http :// archive. nrao . edu/ archi ve/e2earchive .j sp 
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Fig. 2. Correlated flux densities of the source PKS 1402+044 
versus (u, v) distance at 1.6 GHz (top VSOP data), 5 GHz (mid- 
dle VSOP data) and 15 GHz (bottom VLBA data). 



brightness temperature in the source frame in K. The error (lcr) 
is also listed for all the values. The jet shows a wide section be- 
tween 20 and 70 mas (140 - 490 pc projected distance). The 
uniformly-weighted VSOP image shown in Fig. |3f has a higher 
angular resolution (at the expense of considerably higher image 
noise) and indicates that components E and F are essentially re- 
solved. 

The 5-GHz VLBA image in Fig. |3g shows a similar struc- 
ture to the 1.6-GHz VSOP image and the earlier 5-GHz image 
by Gurvits et al. (JT992J. The naturally-weighted VSOP image at 
5 GHz shows that the jet components at 1.6 GHz are resolved 
into a few subcomponents. Here we have differentiated them 
with postfix number in the uniformly-weighted space-ground 
image (Fig. [3]). In this image, the jet appears to be heavily re- 
solved. The core shows a three-component morphology. A weak 
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Fig.3. VLA and VLBI images of PKS 1402+044. The contours are drawn at -2, -1, 1, 2, of the respective first contour (f.c.) 
levels marked in the bottom of each panel. The latter are ~ 3<x of the image thermal noise. The synthesised beams are plotted in the 
bottom-left corner of each image. The basic parameters of each image are listed in Table [2] 



component marked as Al appears at the base of the jet and 
near the brightest component A2. The weakness of component 
Al may be due to synchrotron self-absorption considering its 
high brightness temperature (~ 10 12 K). There are two relatively 
weak jet components, Bl and B2, at 1.6 and 5 GHz between the 
bright components A and C. At the higher frequency, 15 GHz, 
both components are too weak (< 0.9 mJy/beam) to be detected. 
Based on the spectrum at frequencies < 5 GHz, the extrapolated 



total flux density of B1+B2 is ~ 5 mJy at 15 GHz. The nonde- 
tection of the two components indicates that they have a steeper 
spectrum (a < -0.9) at frequencies > 5 GHz. 

The 1.4-GHz VLA image (Fig. [3^) has the lowest resolu- 
tion and shows that there is a weak (~ 33 mJy) component 
(H) at a distance of 3'.'22 and a position angle -107° 4 from the 
core, besides the main emission region. It agrees well with ear- 
lier MERLIN observations made with the Westerbork Synthesis 
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Radio Telescope (WSRT; Gurvits et al. |1992| l. The higher sen- 
sitivity (0.07 mJy/beam) VLA observations (Fig. [3}?) at 5 GHz 
indicate that component H has a weak extension toward east. 
The extension is consistent with the hypothesis that component 
H belongs to the jet of PKS 1402+044. The main emission re- 
gion can be approximated by component G and a combination of 
the inner components (A - F) in the VLA images. Component G 
is also detected at 15 GHz in Fig. [3}; and even 43 GHz in Fig.[3ji. 
The highest observing frequency corresponds to the rest-frame 
emitted frequency of ~ 180 GHz. Arguably, this is one of the 
rare cases of a profound jet emission at millimetre wavelengths. 



4. Discussion 

4. 1 . Spectral properties of the jet 

The resolution (6.88 x 1.19 mas) of the Space VLBI image of 
PKS 1402+044 at 1.6 GHz is close (3.55 x 1.40 mas) to that of 
the ground VLBA image at 5 GHz, enabling extraction of spec- 
tral index information from a combination of the two images. 
We restored the 1.6-GHz VSOP image and the 5-GHz VLBA 
image with a circular Gaussian beam of 4 mas in diameter. The 
artificial beam increases the beam area by a factor of ~ 2 at 
1.6 GHz and ~ 3 at 5 GHz compared to the areas of the original 
synthesised beams. Both images were aligned at the strongest 
component. The shifts in the image centre are less than 0. 1 mas 
(<K 4-mas resolution). After the alignment, the spectral index 
was calculated at all pixels with brightness values higher than 
1.8 mJy/beam (5cr) in the 1.6-GHz image and 1.2 mJy/beam 
(5cr) in the 5-GHz image. A possible core shift between the two 
frequencies as predicted by Kovalev et al. (2008) does not ex- 
ceed 1.5 mas. Thus, it does not affect the large-scale spectral 
distribution. 

The final spectral index distribution between 1.6 GHz and 
5 GHz is displayed in Fig. [4] It shows a smooth distribution of 
spectral index on ~ 20 mas (140 pc) scale. The spectral index 
varies from +0.1 in the optically thick base region to -1.0 in the 
optically thin regions on the western side. To further confirm the 
variation, we plotted the components spectra in Fig. [5] Here we 
also used the 2.3/8.4-GHz visibility data from the USNO RRFID 
database. We fitted the VLBI visibility data with three compo- 
nents at each frequency (1.6, 2.3, 5, 8.4, and 15 GHz). The spec- 
tra of the large-scale components G and H are also plotted. All 
spectra can be approximated by a power-law model, S v = SqV. 
The spectral indices are listed in Table [4] The spectral steepen- 
ing increases with the increase in the distance from the core. The 
spectral difference is 0.47 between components A and B+C and 
reaches 1.57 between component A and the farthest component 

H. This spectral index gradient results in the variation in the flux 
density ratio of the components (B+C) over A from ~ 0.19 at 

I. 6 GHz to ~ 0.05 at 15 GHz. In a sample of sources at various 



Table 2. Parameters of images in Fig. [3] 



Vote 
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Wt. 




f>mi„ 




S peak 


C"rms 


(GHz) 






(mas) 


(mas) 


C) 


(mJy/b) 


(mJy/b) 


1.4 


VLA:A 


NW 


1580 


1310 


7.6 


862 


0.3 


4.8 


VLA: A 


NW 


563 


404 


41.7 


919 


0.07 


15.9 


VLA:A 


NW 


156 


130 


7.8 


754 


0.3 


43.3 


VLA:BC 


UW 


419 


145 


81.2 


476 


0.3 


1.6 


VLBA 


NW 


10.50 


4.85 


1.7 


710 


0.07 




VSOP 


UW 
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1.19 


31.7 


536 


0.3 
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NW 
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1.81 


-4.5 
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0.3 




VSOP 


NW 
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Fig. 4. The spectral index distribution in the jet of the quasar 
PKS 1402+044. 



redshifts, an increase in redshift is equivalent to the increase in 
the intrinsic emitting frequency, v em = v b s (l + z). Thus, a de- 
crease in jet-to-core flux density ratio with increase in redshift 
is expected (Frey et al. [19971 . The quasar PKS 1402+044 has a 
spectral difference of ~ 0.6 at the pc scale, which agrees well 
with the prediction 0.55 + 0.43 by Frey et al. d 19971 1. 
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Fig. 5. The component spectra in the jet of the quasar 
PKS 1402+044. The dotted lines represent the fitted curves with 
the simple power-law model S v — Sov". The estimated spectral 
indices are listed in Table [4] 



4.2. The mass of the central object of PKS 1402+044 

The richness of the core-jet morphology in PKS 1402+044 
makes it a suitable source for estimating parameters of the 
central black hole. The smallest detectable size for a circular 



6 J. Yang et al.: Dual-frequency VSOP imaging of the quasar PKS 1402+044 

Table 3. The parameters and the brightness temperatures of the fitted Gaussian models. 
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Uniformly weighted VLA image at 43.34 GHz 
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Lz 


14 ± 


5 


10.78 ±0.37 


-54.1 ±2.0 


2.21 ±0.73 


0.30 


(3.1 ± 2.3, 


X 10 


Dl 


8± 


3 


10.11 ±0.35 


-73.4 ± 2.0 


1.84 ±0.69 


0.39 


(2.6 ± 2.2) 


x 10 9 


D2 


26 + 


8 


13.18 ±0.38 


-78.0 ± 1.6 


2.52 ± 0.76 


0.22 


(4.3 ± 2.9) 


xlO" 








Naturally weij 


;hted VLBA image at 15.34 GHz 






Al 


102 ± 


13 







0.21 ± 0.02 


0.07 


(7.7 ± 1.8) 


x 10" 


A2 


287 ± 


21 


0.50 ± 0.01 


-19.5 ± 1.3 


0.18 ±0.01 


0.04 


(2.9 ± 0.4) 


x 10 12 


A3 


207 ± 


18 


1.24 ±0.02 


-20.9 ± 0.9 


0.52 ± 0.03 


0.05 


(2.6 ±0.5) 


x 10" 


CI 


29 ± 


9 


9.58 ±0.23 


-46.0 ± 1.4 


1.64 ±0.47 


0.14 


(3.6 ±2.3) 


x 10" 



Table 4. Results of the power-law spectral model fits for each 
component shown in Fig. [5] 



Comp. 


So 


a 


: 1 

X" 




(Jy) 






A 


0.79 ± 0.05 


-0.09 ± 0.04 


1.19 


B+C 


0.17 + 0.02 


-0.56 ±0.10 


0.72 


D 


0.12 ±0.02 


-0.83 ±0.19 


0.16 


G 


0.18 ±0.09 


-0.91 ±0.09 


0.08 


H 


0.06 ± 0.03 


-1.66 ±0.40 





Gaussian component in an image with an rms noise cr rms is de- 
fined as (Lobanov l2005l) : 



<4i 



22-/3/2 



n bmaj b min In 2 In ( 



$ int/^~rms 



t/o-„ 



1/2 



where Si M is the integrated flux density of the component, b m ^ 
and b m i n are the major and minor axes of the restoring beam 
respectively, B = for uniform weighting and/? = 2 for naturally 
weighting. Based on the above criterion, except for the size of 
the combined component from A to F at 43 GHz, all the sizes 
estimated from our VLBI and VLA images in Column (5) of 
Table[3]can be taken as the true sizes of the jet emission regions. 



If the component size is related to the physical transverse 
dimension of the jet, the mass of the central object can be esti- 
mated assuming that the jet is collimated by the ambient mag- 
netic field of the host galaxy. The jet components A2 and A3 
have the best measurements of the width of the jet close to the 
central object, as they are most likely free of the effect of the adi- 
abatic expanding of the jet and synchrotron self-absorption and 
have the higher reliability, djcrd > 10, where d and crj are the 
angular size and error of the fitted circular Gaussian component. 
For a jet collimated by the ambient magnetic field B ext of the 
host galaxy, the mass of the central object Mbh can be related to 
the width of the jet rj et (in pc) according to the following relation 
(Beskin [T997i l: 



M B h - r ]et (B ext /B gI ) 



1/2 



10 13 M R 



(2) 



(1) where B gr is the magnetic field measured at the Schwarzschild 
radius R gr of the central black hole. Equation refers to the 
transverse dimension of the jet measured at distances compara- 
ble to the collimation scale (typically expected to be located at 
10 2 - 10 3 R gl ). A typical galactic magnetic field is B e xt ~ 1CT 5 G 

10 4 G (Field & 



(Beck 2000) and one can expect to have B gr 
Rogers [T993l) . Based on these parameters, the mass of the central 
object is ~ 1O 9 M . The main uncertainty of the mass estimation 
arises from the uncertainty in B gr . However, the dependence of 
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the mass on the value B gr is rather weak, Mbh <* B g , • , with the 
magnetic field varying within 4 orders of magnitude, the esti- 
mated central black hole mass varies within two orders of mag- 
nitude. The external magnetic field Z? ext normally varies within a 
narrow range around (10~ 5±I G). Thus, the magnetic field uncer- 
tainty should not affect the estimated mass drastically. 



4.3. Brightness temperature 

Based on the parameters of the Gaussian models listed in 
Table [2] we calculated the brightness temperature of each com- 
ponent using the following formula (Kellermann & Owen 1988 ): 



1.22xl0 12 (l+z)4^ 



(3) 



where S mt is the integrated flux density in Jy, d the size of a 
circular Gaussian component in mas, and v the observing fre- 
quency in GHz. The estimated brightness temperatures are listed 
in the last column of Table [2] Among these components, the 
component A2 has the highest brightness temperature, T B = 
(5.7 ± 0.5) x 10 12 K, which is somewhat higher than the inverse 
Compton limit (~ 10 12 K; Kellermann & Pauliny-Toth [19691 
but still 10-times lower than the currently known highest value 
of 5.8 X 10 B K f ound in the BL Lac object AO 0235+164 by 
Frey et al. (120001) . In the equipartition jet model of Blandford 
& Konigl j 19791 1, the limiting brightness temperature is about 
3 x 10 n <5 5 / 6 K, where 5 is the Doppler factor. Comparing the 
theoretical value with the estimated brightness temperature of 
the component A2, we can obtain a conservative lower limit to 
the Doppler factor of the inner jet, 5 * 23.7. 

The variation in the observed brightness temperature with 
increasing distance from the core is plotted in Fig. [4] Following 
Marscher ( 1990), we assume that each of the jet components is 
an independent plane shock in which the radio emission is dom- 
inated by adiabatic energy losses. The jet plasma has a power- 
law energy distribution, N(E)dE oc E~ s dE. The magnetic field 
varies as B oc where / is the distance from the central object. 
The Doppler factor is assumed to vary weakly throughout the jet. 
Under these assumptions, one can relate the brightness tempera- 
ture Tbjet of each jet component to the brightness temperature of 
the core r b , C ore: 



^b,jet — ^b,core(^jet/^core) 



(4) 



where d represents the measured size of the core and jet features 
and e = [2(2s + 1) + 3a(s + l)]/6 (Lobanov et al. |200TT ). We 
take s = 2.06 corresponding to the synchrotron emission with 
the spectral index of component B+C a = 0.53, and a — 1 corre- 
sponding to the transverse orientation of magnetic field in the jet 
(Lobanov et al. 2001). At each frequency, we take the brightest 
component as the core. Comparing with the measured one, we 
plotted the ratio of r o bs/T t h e ory versus the distance. The largest 
discrepancies occur with components B at 1.6 GHz, Bl, and B2 
at 5 GHz. These components have ratios d/crj = 5 at 1.6 GHz 
and djcTd < 3 at 5 GHz. The ratios indicate that the discrep- 
ancy may be caused by uncertainties in the size estimates for 
strongly resolved components or inhomogeneities in the plasma 
over such large emitting regions. 

4.4. Proper motion 

Using an earlier VLBI observation at 5 GHz in 1986 by 
Gurvits et al. (1992), we tried to estimate the proper motion in 
PKS 1402+044 over the time interval of 15 years. We also fitted 



1E14 
1E13 
1E12 
1E11 
1E10 
1E9 
1E8 



A 

f 



-•4- 1.6 GHz 
-<S>-- 5 GHz 



B 



C 

--■I D 

ff" 



I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I I I I I | I I I I I 1 I I 

10 20 30 40 50 



O 



E 

"A 



4-Ll. 



10 20 30 40 50 60 

Projected distance (mas) 

Fig. 6. Brightness temperature variations along the jet in PKS 
1402+044 (top) and comparison with the model predicted value 
(bottom, Marscher 1990). The component identification is de- 
scribed in Fig. [3] The triangles represent the estimated bright- 
ness temperature at 1.6 GHz; the circles represent the brightness 
temperature at 5 GHz. 



the (m,v) visibility data with 6 Gaussian models and found that 
the position offset is within the one fifth of the beam (10 x 2 
mas in PA. -2°) of the image in 1986 in the east-west direction 
where the image has the higher resolution, and there is no con- 
sistent shift. Thus, an upper limit of the apparent proper motion 
li in the EW direction is 0.03 mas yr _1 . This corresponds to the 
apparent velocity upper limit of B- dpp = 3c based on the relation 
/J app = 1.58 x 10~V) A (1 + z) (Kellermann et al. 120041 , where 
the angular size distance Da is measured in Mpc, [i in mas yr~' 
and y8 ap p in the unit of the speed of light, c. 

Using the determined lower limit S = 23.7 (6 » sJbI vp + 1) 
and the following equations (e.g. Hong et al. 120081 ): 



J = 



tan® = 



lapp + S 1 + 1 

25 

2/?app 
£app + <5 2 -l' 



(5) 
(6) 



a lower limit to the Lorentz factor y » 12 of the jet and an upper 
limit to the viewing angle to the line of sight <p ~ 1° can be 
determined. All the estimates suggest that PKS 1402+044 is a 
relativistically beamed radio sources. 



5. Summary 

Based on multi-frequency VLBI (1.6, 2.3, 5, 8.4, and 15 GHz) 
observations including dual-frequency (1.6 and 5 GHz) VSOP 



8 



J. Yang et al.: Dual-frequency VSOP imaging of the quasar PKS 1402+044 



observations and VLA (1.4, 5, 15, and 43 GHz) observations of 
the high-redshift quasar PKS 1402+044, we draw the following 
conclusions. 



well-defined 
-23 kpc from 



1. The quasar PKS 1402+044 demonstrates a 
core-jet morphology that can be traced out to 
the source core. 

2. The radio spectral index distribution and the component 
spectra prove that the jet has the steeper spectrum with in- 
creasing distance from the core. 

3. Based on the measurement of the transverse size of the jet, 
and assuming that the external magnetic field collimates the 
jet model, the mass of the central object is estimated as ~ 
1O 9 M . 

4. PKS 1402+044 has a bright core (5.7 x 10 12 K), and the ob- 
served brightness temperature variation is basically consis- 
tent with the shock-in-jet model. 

5. No firm detection of a proper motion in the jet can be made. 
An upper limit of the apparent proper motion in the east- 
west direction is 0.03 mas yr _1 , corresponding to the appar- 
ent speed of 3 c. 

Based on the lower limit of the Doppler factor 6 = 23.7 (6 » 
'ySapp + 1), we estimate the lower limit to the Lorentz factor 

y=12 and the upper limit to the viewing angle of the inner 
jet to the line of sight as <p = 1°. 
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